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ABSTRACT 

The present document is a comprehensive and fairly self-contained study 
of centrally obscured optical transmitting and receiving antennas and is intended 
for use by the laser radar- and communication systems designer. The material 
is presented in a format which allows the rapid and accurate evaluation of antenna 
gain. The i resnel approximation to scalar wave theory is reviewed and the 
antenna analysis proceeds in terms of the power gain. Conventional range equa- 
tions may then be used to calculate the power budget. 

The transmitter calculations, resulting in near and far field antenna gain 
patterns, assumes the antenna is illuminated by a laser operating in the fundamental 
cavity mode TEM 00 . A simple equation is derived for matching the incident 
source distribution to a general antenna configuration for maximum on-axis gain. 

An interpretation of the resultant gam curves allows a number of auxiliary design 
curves to be drawn which display the losses in antenna gain due to pointing errors 
and the cone angle of the outgoing beam as a function of antenna size and central 
obscuration. The use of telescope defocusing as an approach to spreading the 
beam for target acquisition is compared to some alternate methods. 

Incorporated into the document is an analysis of optical receivers which 
includes a discussion of the optical losses at the detector. Direct and heterodyne 
detection schemes are considered separately. Included is a detailed discussion 
on the manner in which the local oscillator field distribution across the detector 
affects the overall receiver gain in heterodyne systems. 
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TRANSMITTER AND RECEIVER ANTENNA GAIN 
ANALYSIS FOR LASER RADAR 
AND COMMUNICATION SYSTEMS 


1. INTRODUCTION 

Laser communication and tracking systems have long been predicted to have 
immense capacity to supplement and extend the capabilities of the more conventional 
microwave and mill' meter wave systems. These projections are now being 
realized in the form of laser tracking systems which have made great advances 
in the accuracy of orbit determination for both man-made earth satellites and 
the moon, as well as in the measurement of point-to-point terrestrial distances. 
Experimental laser communications systems have been developed and systems 
studies indicate that lasers may represent the optimum approach to transferring 
high data rates from earth observational spacecraft to NASA's planned Tracking 
and Data Relay Satellite System (TDRSS). All of these applications share a 
common analytical base, the range equation for either radar or communications. 

One of the most important parameters in the range equation is the antenna gain, 
for both transmitting and receiving antennas. This document had its origin in 
the recognition that the data available on optical transmitting and receiving 
antennas was too sparse and incomplete to provide a reasonable and sufficiently 
accurate basis for system design. In addition, much of the available information 
was widely scattered among the various scientific and engineering journals and 
often presented in a format not directly useful to the systems engineer. It 
became obvious that it would be most desirable to generate a set of comprehensive 
systems designer-oriented plots which would permit the rapid and convenient 
determination of transmitting and receiving antenna gain. That has been the authors’ 
objective in ■writing this document. 

Laser communications and radar systems employ conventional telescopes 
as optical antennas. These telescopes frequently have central obscurations, such 
as the secondary mirror in the Cassegrainian telescope shown in Fig. 1. It was, 
therefore, necessary to examine the effects of such obscurations. The narrow 
beams (i.e., high antenna gains) associated with laser communications and radar 
systems require precise pointing and the systems designer must be able to 
relate pointing errors to losses in antenna gain. In addition, the narrow beams, 
when coupled with the finite velocity of light, necessitate a point-ahead angle 
which must be considered by the systems designer; this imposes a further require- 
ment that the systems engineer know the distribution of antenna gain with off-axis 
angle. A knowledge of the antenna gain distribution is of more immediate usefulness 
than the normalized intensity plots previously recorded by Campbell and DeShazei 1 
for the near field and by Buck 2 and Kaufman 3 for the far field. 
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Figure 1. Cossfcgrainian Telescope 


The establishment of communications or tracking in a laser system requires 
that an acquisition sequence be performed, again because of the narrow beams 
which are usually much narrower than the angular uncertainties which must be 
coped with by the system. One approach to acquisition is to utilize a diverged 
transmitter beam, much wider than the ultimate beamwirttb required by bit 
error rate or tracking accuracy considerations, which is diverged to the maximum 
extent possible within the limits of the required acquisition signal-to-noise ratio. 
One means to diverge the transmitter Learn is to defocus the optical antenna 
which effectively transforms the far field pattern to a near field distribution. 
Defocusing by arbitrary amounts can, however, result in a local intensity 
minimum on-axis in the i'ar field which may be undesirable for acquisition. 

Thus, the effect of defocusing must be considered in detail. 

In the present document, we review the fundamental assumptions inherent 
in the Fresnel approximation to scalar wave theory and apply the theory to 
centrally obscured circular apertures illuminated by Gaussian amplitude wave- 
fronts of varying curvature using a model described in Fig. 1. The beam is 
characterized by a beam width w measured from the axis to the l/e 2 intensity 
point and a phase front of curvature R. Both quantities are measured in the plane 
containing the secondary mirror. The telescope system is described by two 
parameters, a primary radius "a" and a secondary radius "b". The on-axis 
power gain is then computed as a function of the various parameters in both the 
near and far field, and an expression is derived which yields the optimum aperture 
radius to beam width ratio (a/o>) for maximum on-axis gain in the far field as 
a function of obscuration ratio (b/a). The off-axis near and far field power gain 
distributions are then computed for the optimum antenna configurations. A 
series of auxiliary design curves allow the rapid evaluation of the on or off-axis 
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gain, the truncation loss by the aperture, and the effective cone angle for an 
arbitrary wavelength and aperture size. A comparison is made between unapertured 
laser beams and defocussed antennas when large cone angles are required for 
acquisition. Finally, the use of centrally obscured antennas as receivers is 
considered. Calculations of the receiver gain are made which include computa- 
tions of the losses due to blockage of the light by the central obscuration and 
its effect on the field distribution at the detector. The losses incurred by systems 
utilizing direct detection and heterodyne (or homodyne) detection are separately 
discussed and the procedure for matching the detector dimensions to the telescope 
system for maximum efficiency is considered in detail. 


2. THE FRESNEL DIFFRACTION FORMULATION 

In the present analysis, we choose spherical coordinates for the calculation 
of the field distribution at the point of observation (r ^ , 6 , </> ). Figure 2 shows 
the diffracting aperture to be in the $ 0 = tt/2 plane and a point in the aperture 
designated by the coordinates ($, ,tt/2, <£„). 



Figure 2. Coordinate System 


The distance between points, r 01 , also shown in Fig.2, is obtained from 


r oi = ( x i “ V 2 + <Vi - ^o ) 1 + (*i ) 2 

where z 0 = 0. In spherical coordinates, the latter expression becomes 


r oi " ( r | Sin 0 t Cos0j - r 0 Cos<^ 0 ) 2 + (r t Sintfj Sin^ - r 0 Sin<£ 0 ) 2 + (r l Cos$ l ) 2 

' "i.' 
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which reduces to 


r oi = r i 


1 + _ - 2!2.Sin0j Cos^ - 4 > 0 ) 


1/2 


The Fresnel approximation assumes that the first two terms of the binomial 
expansion of the square root are adequate. Thus we have 


r 


01 



— - Sin6>j Cos(^>j - < i> ) 

2r ? l 


( 1 ) 


Huygen's principle states that each point in the aperture is the source of a 
spherical wave weighted by the amplitude of the field at the point. The electric 
field at the observation point (r x , 9 ^ , 4>^ ) is thus given by an integration over 
the aperture, that is, 4 


E ( r i> e v *i> = J* J e)krn ' E o( r o- d ^o r o dr o 


( 2 ) 


where a is the aperture radius, X is the wavelength, and k is the propagation 
constant (277 A). If the field distribution at the aperture is circularly symmetric, 
then Eq. (2) may be written 


E ( r i • e v 4> x ) - 


A r i 


•i r a ci* 
: J 1 


jkr 


0 /2r ! 


E n( r o) 


•jk Sinfij Cos(0j"^ o ) 


d ^o r o dr o 


( 3 ) 


where we have used Eq. (1) to approximate r Q1 in the exponential argument in 
(2) and the approximation r 0 j ~ r t in the less critical denominator term. The 
circular symmetry of our system requires the field at the observation point to 
be independent of <p l > and we may, therefore, choose a value of 4 > l for convenience. 
Choosing </>, = - n/2 and using the integral representation of the Bessel function 
of the first kind given by 5 
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we may write Eq. (3) as 


i r 2n 

1 p i ( z Sin <£ — n <t>) 

277 1 
•'O 


E(r r 0,) = 


-jk jkr, f 

” J„ 


l ik, 0 '' Jr, J 0 (l<r 0 Sin(J 1 )E 0 <r 0 ) r„dr„. 


The latter equation is a valid approximation if the condition 


(4) 


j£(- tJ |S 1 n # 1 iy < <l 

is met. The above expression is derived from the third term in the binomial 
expansion for r Q1 and sets limits on the observation regions for which the 
Fresnel approximation is valid. 


3. TRANSMITTER GAIN 

3.1 Definition of the Optical Antenna Gain 

We assume the incident source function E 0 (r 0 ) to have a Gaussian amplitude 
function (corresponding to the fundamental laser mode) and tc be given by the 
equation 


V r o>-- 



e 


jkr 7 /2R 


where o is the distance from the axis to the l^e* intensity point and R is the 
curvature of the phase front at the aperture (See Appendix A). The proportionality 
constant normalizes the incident power P over the entire aperture plane to one, i.e., 


P ; 




V lr o d ^o 


1 . 


(5) 
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We assume this field function E Q (r Q ) to impinge upon a centrally obscured ap- 
erture having an opaque obscuring element of radius b. The impinging field is 
centered about the telescope axis. Axicon devices and off-axis incident fields 
will not be considered in the present work. The latter approaches attempt to 
redistribute the field in order to minimize the obscuration by the secondary 
mirror and the resultant degradation in transmitted power. Peters and Ledger 6 
consider axicon devices in great detail and note the advantage in particular for 
an obscuration ratio (b/a) of 0.5. This central obscuration to aperture ratio is 
larger than is commonly employed. For smaller b/a values, the improvement 
in on-axis gain is marginal and does not appear to justify the more stringent 
alignment tolerances. 7 

With the substitution of E 0 (r 0 ) in Eq. (4), the intensity 8 distribution at the 
observation point (r } , 9 ^ ) is given by 


I(r lt Q ) = 


k 2 



/2 

■/tTCx) 


-'0 2/ « 


1 i 
e 



2 



2 


J 0 ( kr Q Sin 8 y ) r Q dr 0 


( 6 ) 


The intensity of the radiation field produced at a distance rj from a unit power 
isotropic radiator is, of course, given by 


Io = 


4vr 2 


In the ensuing calculations, we choose to define the gain by the expression 


I(r., 0 .) 

G(v e t) = — - — - 


~ |j e 2 ^ ^ e roA> Jo( kr o Sin 0 i) r o dr o 


(?) 


where the gain depends on r t only in the near field. 
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It is important to note that the above definition of gain differs from the 
definition conventionally used in microwave theory 5 in that the power of the 
. isotropic radiator is set equal to the total power (unity) impinging on the aperture 

T plane as opposed to the power actually transmitted by the aperture. The latter 

is given by 


770 )' 


rr 

/ n A 


2 2 
-2 rt/u 


r o dr o d ^o 
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- 2 b 2 /<» 2 _ e -2 f . 2 /*) 2 

= e e 


< 1. 


Thus, the gain calculated under the present definition will be smaller by the 
factor (e" 2 “ 7 - e~ 2a ) where a= aA> and y = b/a. Our choice of definition 

is motivated by the fact that, for the systems presently under discussion, 
truncation by the aperture boundary and obscuration by the secondary represent 
real losses in transmitter power and hence should be included in a calculation 
of the overall antenna gain. From the conventional definition of gain, we 
would obtain the result that a !<x> - 0 (uniform illumination) is optimum for maxi- 
mum on-axis gain in the far field. This is a realistic result only if all of the 
available energy is actually transmitted by the aperture. Therefore, the con- 
ventional definition might be appropriate for systems where the energy is re- 
distributed by means of an axicon or other device; it is not appropriate for the 
systems presently under discussion. 

The qualitative behavior of the curves in Fig. 3 can be explained quite 
easily. As a very large spot sizew is made smaller relative to the aperture 
radius (aA» increasing), less energy is lost due to truncation and the fractional 
power transmitted increases. For centrally obscured antennas (y ^ 0), however, 
a further reduction in the spot size results in a significant fraction being blocked 
by the secondary with the net result that the curves corresponding to obscured 
apertures exhibit a minimum loss point. As we shall see later, however, the 
a/«- ratios corresponding to maximum fractional transmission by the aperture 
do not coincide with those values which give the maximum on-a>is gain. 

Returning to Eq. (7) we make a substitution of r 0 - s«a, where s is a dimen- 
sionless variable and a is the aperture radius, and obtain 
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Figure 3. Combined trunculion and obscuration loss in aB 
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We now define a number of dimensionless parameters to simplify L'q. ^s) and 
rewrite it as 


G x (u. » , X) - 8 



a 2 a 2 


f ’ ei>*“ . 

S 2 


1 


''** 'M 0 (X u) - 2 du 


(!» 
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where we have defined 


u 


s 


2 


C? / Cu 


y - b/a 
X = k a S i n 9. 


ka 2 /l 1 


2 \r, R 


(10) 


Eq. 9 can be rewritten us 


V a - 5, 7 , X) - R ,(a, A >. X) 

\2 1 


( 11 ) 


where u and , are system parameters X - ka sm is a convenient plotting 
parameter for off-axis distributions, - c includes near field and defocusing effects 
aud A - •• a 2 is the aperture area. The factor <1 A A 2 is the well known upper 
limit on antenna gain which is obtained for a uniformly illuminated unobscured 
circular aperture." 


In the following sections, we examine the function g(.-< , X) which we 
shall refer to as the "transmitter efficiency factor" and which is defined by 
the equation 


k. 

m * 


. 7 . X) 



o~ a "LfX- u) du‘ 


( 12 ) 


2.2 Un-Axis Gain 

A general equation for the on-axis antenna efficiency for incident Gaussian 
energy may be obtained by setting X 0 in Eq. (12) and performing the rather 
simple integration to obtain 






3C 



r>* 


X 


r 2 a- + e -*«V _ 2e - a W + ^ cos L3(y 2 - 1)]} (13) 

If we specialize to the far field and Gaussian amplitude plane waves (.5 = 0) 
impinging on an unobscured aperture (y = 0), Eq. (13) becomes 


g T (u, 0, 0, 0) = — [e-* 2 - 1] 2 . ( 14 ) 

a 2 

Setting the derivative of Eq. (14) equal to zero yields the transcendental 
equation 



[2cg * fl-e"** = 1 ( 15 ) 

which can be solved numerically or graphically to obtain = 1.12 for the 
aperture to spot size ratio that yields the maximum on-axis gain for the un- 
obscured case. The corresponding value of the efficiency factor is equal to 
0.8145. 


For a non-zero obscuration ratio, Eq. (13) becomes 


0, y, 0) 




(16) 


If we take the derivative with respect to a , we again obtain a transcendental 
equation, that is 


\ 


t 


4* 



2a 2 + 1 
2a 2 y 2 + 1 


e -a (1 


■r 2 ) = 


1 


(17) 


An approximate solution obtained using second order perturbation theory (See 
Appendix B) is given by 


i 


a » 1 . 12 - 1 . 30 >' 2 + 2 . 1 2 y 4 ■ ( 1S ) 

Eq. (18) gives the optimum aperture to beam width ratio for a general obscuration 
aud is accurate to within ±1% for y < 0.4. 
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Figure 4 illustrates the degradation in the far field efficiency factor experi- 
enced when non-optimum antenna configurations are used. Each curve is labeled 
by a particular value of obscuration ratio y , and the efficiency factor is plotted 
(in dB) versus the aperture radius to beam width ratio. 





Figure 4. Far field axial gain relative to 4/rA/)? as a function of a. for 

five different obscuration ratios ( y ) I 

|. 

W- 

We now consider the on-axis gain in the near field for the optimum antennas 

just derived. This can be obtained by choosing a particular value for y , calculating 

the corresponding value of a from Eq. (18), substituting both values into Eq. (13) 
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and plotting the efficiency factor versus ,8. Figure 5, where we have chosen to 
plot 10 Log g (a . , 8 , 7 , 0) versus ,6 / 2 tt, clearly displays a series of gain minima. 

As we shall see in the next section, these minima correspond to observation 
regions where a local intensity minimum (shadow) occurs on-axis in the gain 
profile. It is worthwhile to note that the minima occur for unobscured as well 
as obscured antennas and, therefore, are not related to the shadow produced by 
the secondary mirror. Instead they result from destructive interference between 
the Huygens wavelets. Their positions correspond approximately to the well-known 
intensity nulls which are observed at discrete distances (integer values of S/2n) 
from a uniformly illuminated aperture . 1 In the case of uniform illumination, 
the cancellation of electric field is total; Gaussian illumination produces only 
partial cancellation. In addition, the positions of the minima are displaced 
slightly in going from a uniform to a Gaussian distribution. Increasing the 
obscuration ratio from zero displaces the minima toward larger values of 3. 

For plane waves incident on the aperture, the far field value of the on-axis 
efficiency factor is, of course, the value at /3 - 0. As the distance between the 
aperture and the observation plane decreases, ~ increases. From our definition 
for however, other information is obtainable from Fig. 5. The effect of 
phase front curvature at the aperture on the efficiency factor in either the near 
or far field can be determined by calculating the appropriate value for 8 (includ- 
ing both near field and phase curvature terms as given by equation ( 10 )) and 
reading the corresponding efficiency factor from Fig. 5. Thus, Fig. 5 is a 
concise method for describing an infinite number of experimental situations. 


3.3 Antenna Gain Patterns 

We now turn our attention to the off-axis gain {0 ^ 0). Figures 6 through 10 
are computer generated plots of the transmitter efficiency factor given by Eq. 
(12). Details pertaining to the numerical techniques employed in the calculation 
are included in Appendix C as is a complete listing of the computer program. 
Each figure assumes specific values for 7 and a. The latter parameter is cal- 
culated from Eq. (18) for maximum axial gain in the far field. The curves dis- 
played within each figure are labeled by the appropriate values of fi/2v. In part 
(a) of Figures 6 through 10 values of j3/2rrwere chosen corresponding to values 
midway between adjacent extrema (maxima and minima) in Fig. 5. In part (b), 
values corresponding to the first few extrema in Fig. 5 are considered. Com- 
paring Figs. 6 through 10 with Fig. 5 verifies the assertion made in the previous 
section that the minima in Fig. 5 correspond to local intensity minima on-axis 
("shadows" in the central lobe). The maxima in Fig. 5, on the other hand, are 
associated with localized enhancement of the gain near the axis. 
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Figure 5. Axial gain relative to 477A/>? as one proceeds from the far field = 0) to the near 
field (/}/ 0) for four different obscuration ratios. The same curve can be used to determine the 
degradation in the far field gain due to defocusing (see e q. (10) in the text). The minima core- 
spond to shadows in the central lobe of the antenna pattern. These shadows occur at certain dis- 
tances from the aperture in the near field and also in the far field when the system is sufficiently 
defocused. 



10LOG[g T («, fi, y, X)] 


Ml] 


-20.0-16.0-12.0 -8.0 -4.0 0.0 4.0 

X=kaSir>0 t 


12.0 16.0 20.0 


Figure 6a. Far field {/3» 0) and near field 0) transmitting antenna gain distributions in dB 
relative to 47 tA/X? as a function of the angle 6, from the optical axis of the antenna. The aper- 
ture to gaussian spot size ratio a is chosen to give maximum on^xis gain in the far field accord- 
ing to Eq. (18) in the text. The near field plots are labeled by the value /?'= /3/2nand also cor- 
respond to fcr field gain distributions obtained from defocused antennas. 
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Figure 6b. Near field (or defocused far field) transmitting antenna gain distributions in dB rela- 
tive to 4rrA/}? as a function of the angle $, from the optical axis of the antenna. The aperture 
to gaussian spot size ratio ais chosen to give maximum on-oxis gain in the far field accor ,g to 
Eq, (18) in the text. The near field plots are labeled by the values /3' * fiflTT correspon j ; ng to 
the maxima and minima in Fig. (5). Comparisons with the latter figure indicate that the minima 
in Fig. (5) correspond to on-axis shadows in the central lobe and that the maxima correspond to 
localized enhancements of the on-axis gain. 


15 



-5.0 


/ J '=0 
/?' = 0.53 
fl' = 1.23 


—15.0 


- 20.0 


-25.0 


-30.0 


-35.0 


-40.0 

-20.0-1.6.0-12.0 -8.0 -4.0 0.0 4.0 8.0 12.0 16.0 20.0 

X = kaSin0, 


Figure 7o. For field (/3 = 0) and near field (/3/ 0) transmitting antenna gain distributions in dB 
relative to 47 tA /\ 2 as a function of the angle 6, from the optical axis of the antenna. The aper- 
ture to gaussian spot size ratio <x is chosen to give maximum on*oxis gain in the far field accord- 
ing to Eq. (18) in the text. The near field plots are labeled by the value /?' = ^/2rrand also cor- 
respond to far field gain distributions obtained from defocused antennas. 
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Figure 7b. Near field (or defocused far field) transmitting ontenna gain distributions in dB rela- 
tive to 4<7A/\^ cs a function of the angle Q, from the optical axis of the antenna. The aperture 
to gaussian spot size ratio ais chosen to give maximum on-axis gain in the far field according to 
tq. (18) in rhe text. The near field plots are labeled by the volues £' = {$/7w corresponding to 
the maxima and minima in Fig. (5). Comparisons with the latter figure indicate that the minima 
in Fig. (5) correspond to on-axis shadows in the central lobe and that the maxima correspond to 
localized enhancements of the on-axis gain. 
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Figure 8a, Far field (/ 3 = 0) and near field (/3W 0) transmitting antenna gain distributes in dB 
relative to 4 ttA/A? as a function of the angle 6 t from the optical axis of the antenna. The aper- 
ture to gauss ian spot size ratio a is chosen to give maximum on-axis gain in the far field accord- 
ing to Eq. (18) in the text. The near field plots are labeled by the value /3* = (Y2ti and also cor- 
respond to far field gain distributions obtained from defocused antennas. 
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Figure 8b. Near field (or defocused far field) transmitting antenna gain distributions in dB rela- 
tive to 4 77 A/ X? as a function of the ^nqie 6, from the optical axis of the antennc. The aperture 
to gaussian spot size ratio ais chosen to give maximum oncxis gain in the far field according to 
Eq. (18) in the text. The near field plots ore labeled by the values / 3 # = /3/2tt corresponding to 
the maxima and minima Fig. (5). Comparisons with the latter figure indicate that the minima 
in Fig. (5) correspond \o on-axis shadows in the central lobe and that the maxima correspond to 
localized enhancements of the on-oxis gain. 
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Figure 9a. Far field (.6’= 0) and near field (.3/ 0) transmitting antenna gain distributions in JB 
ielative to 4 t 7A/A? as a function of the angle 0, from the optical axis of the ontenna. The aper- 
ture to gaussian spot size ratio a is chosen to give maximum on-axis gain in the far field accord- 
ing to Eq* (18) in the text. The near field plots are labeled by the value fi e = ,o/2~ 7 and also cor- 
respond to far field gain distributions obtained from defocusea antennas. 
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Figure 9b. Near field (or defocusud far field) transmitting antenna gain distributions in d8 rela- 
tive to 4' r ?A/ / \^ as a function of the angle ** # from the optical axis of the antenna. The aperture 
to gaussian spot size ratio a is chosen to give maximum on-axis gain in the far {.eld according to 
Eq. (18) in the text. The near field plots are labeled by the values ' -V *2 ^ corresponding to 
tne maxima and minima in Fig. (5). Comparisons with the latter figure ind : cate that the minima 
in Fig. (5) correspond to on-oxis shadows in the central lobe and that the maxima correspond to 
localized enhancements of the on-axis gain. 
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Figure 10a. Far field (,~ = 0) and near field (/V 0) transmitting antenna gain distributions in dB 
relative ro 4 77 A/ A 2 as a function of the angle 0, from the optical axis of the antenna. The aper* 
ture to gcjssian spot sue ratio a is chosen to give maximum on-axis gain in the far field accord** 
ing to Eq. (18) .n the text. The noar field plots are labeled by the value /3* - /Y2ti end also cor* 
respond to far field gain distributions obtained from defocused antennas. 
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Figure 10b* Near field (or defocused far field) transmitting antf/ina gain d isfributions m dB relc* 
live to 4' r /A/A? as a function of the angle 6 t from the optical axis of the antenna. The aperture 
to gaussian spot size ratio ais chosen to give maximum on-axis gain in the far field according to 
Eq. (18) in the text* The near field plots are labeled by the values /V = fvlv corresponding to 
the maxima and minima in Fig. (5)* Comparisons with the latter figure indicate that the minima 
in Fig. (5) correspond to on-axis shadows in the central lobe and that the maxima correspond to 
localized enhancements of the on-axis gain. 



In all the curves, the efficiency factor is plotted against the dimensionless 
parameter X = ka Sin f,. This is a particularly useful format since the spreading 
of the bean* in the far field and the dependence of the cone angle on aperture 
size is automatically taken into account. This eliminates cumbersome scale 
transformations in going from one observation distance to another or to different 
aperture sizes or wavelengths. 


3.4 Transmitter Summary 

It is now possible to generate several useful curves which greatly simplify 
the evaluation of optical antennas for a particular application. From Eq. (11), 
the transmitter gain in dB for a general antenna at an arbitrary observation 
point is given by 


G T (db) = 10* Log 


The first term of Eq. (19) is plotted in Fig. 11 for several important wavelengths 
as a function of aperture radius. The second term, corresponding to the efficiency 
factor in dB, is a function of four variables. The on-axis (X = 0) efficiency factor 
for optimized focused systems ( ,i given by Fq. (18)) in the far field (.: = 0) is 
plotted versus in Fig. 12. Near field or dei'ocusing effects ( ' / 0) on the 
efficiency factor for these optimized antennas can be estimated reasonably well 
using Fig. 5. Far field efficiency factors for arbitrary non-optimum antenna 
structures can be obtained from Fig. 4. 

The off-axis (X / 0) gain normalized to the on-axis gain is displayed in 
Fig. 13 for four optimized antennas corresponding to ■ = 0.0, 0.2, 0.3, and 0.4. 
These curves indicate the dB loss in far field gain attributed to pointing error 
and are plotted versus the angular error in units of \ /D radians where D is 
the antenna aperture diameter. Fig. 13 can in turn be used to define the cone^ 
angle of the outgoing beam. If we choose to define the half angle of the cone, J , 
as the angle from the optical axis to the 1/e 2 point on the power gain distribution 
(Fig. 13), we obtain the curves presented in Fig. 14. In the latter figure, each 
shaded region is associated with its corresponding laser wavelength. For a given 
aperture radius, the upper and lower limits of the shaded regions are defined 
(somewhat arbitrarily) by the angular spread of the optimized antennas for the 
unobscured case / (y = 0.0) = 0.94 \ /Dj and for an obscuration ratio of 0.4 
= 0.4) = 0.79 \ /Dj respective’y where the latter relations were obtained 
from ihe -8.7 dB (1/e 2 ) points in Fig. 13. One can interpolate between these 
two limits to obtain cone half-angles for intermediate values of y . 
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Figure 14. The half angle of the far-field transmitter cone produced by an antenna configuration 
optimized to give maximum on-saxis gain as a function of the aperture radius a and wavelength K 
The half angle is defined as the angle from the optical axis to the 1/e 2 intensity points. 




Clearly, the total gain of an optimized antenna at any point in the far field 
can be obtained by simply summing the separate contributions from Figs. 11, 

12, and 13, or from Figs. 6 through 11. The near field gain for optimized antennas 
is obtained by summing the contributions from Figs. 5 and 11 or by using Figs. 6 
through 11 depending on whether the observation point is on or off the optical axis. 
The on-axis far field gains for arbitrary non-optimum antennas are calculated 
by summing the contributions from Figs. 4 and 11. For systems where the 
conventional definition of gain as described in Section 3.1 is more appropriate 
for analysis, an additional corrective factor is obtained from Fig. 3 and added 
to the total gain. The latter figure also gives the dB loss in total transmitted 
power due to the combined effects of truncation and central obscuration for 
arbitrary a and y . The individual losses due to the aperture boundary and 
obscuration can be obtained by reading the truncation loss from the curve 
corresponding to y = 0 in Fig. 3 and subtracting that value from the combined 
loss. 

Figures 6 through 10 indicate that defocusing can indeed be used to spread 
the beam continuously over larger angles in tb inr ..eld. Generally as the 
obscuration ratio increases, the gr\. profile tends to become more undulatory 
with angle. This lack of "smoothness" in the profile does not prevent the use 
of a defocused antenna in the acquisition of a target, however, since most acquisi- 
tion techniques merely require the presence of a "star" on which to lock the 
tracking servos. For discrete jumps to cone angles on order of magnitude larger 
or more, one might consider transferring the laser beam to a mucn smaller 
parallel antenna or bypassing the antenna optics altogether and utilizing the 
natural divergence of the laser beam. The on-axis gain in dB of an untruncated 
laser beam is given by 


G(db) = 10- Log 



( 20 ) 


where % is the radius of the laser beam waist as defined by Siegman 10 and is 
typically on the order of a few millimeters. Eq. (20) is plotted as a function of 
waist radius for several notable wavelengths in Fig. 15. The beam divergence 
in the far field for an untruncated laser beam is given by the expression 10 



(21) 
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and is plotted in Fig. 16. The use of smaller parallel antennas or untruncated 
laser beams would result in much smoother gain profiles in the far field ~ 
especially when compared to highly defocused (fi = 100) obscured antennas which 
exhibit large excursions in the far field intensity of 10 dB and more near the 
optical axis as in Fig. 17(b). For systems operating on marginal signal to noise 
ratios, this would be equivalent to "blind" spots in the overall antenna pattern. 

In concluding the transmitter discussion it is worthwhile to mention that a 
general computer program has been developed by R. F. Schmidt 14 which utilizes 
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vector diffraction theory and is capable of computing the effect of misalignment and 
secondary mirror supporting struts on the Fresnel and Fraunhofer patterns. The 
program has been used successfully in the evaluation of several microwave 
systems but would require modifications for application to optical and near 
infrared wavelength systems dve principally to the integration algorithms used. 



Figure 17a. Far field gain relative to 47TA/A. 2 for a highly defocused antenna (/3^= 100) when the 

antenna is unobscured. 
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3.5 An Illustrative Example 

Consider a Cassegrainian optical antenna having a 10 cm exit aperture 
radius a, a 2 cm secondary mirror radius b, and illuminated by a laser operating 
in the Gaussian mode at 10.6 microns. The value of y corresponding to this 
antenna is thus y =0.2. The far field case with plane waves incident on the 
antenna (,6 = 0) will be illustrated first, where we assume that the maximum on- 
axis gain from the antenna would be desirable. A value for a of 1.07 is obtained 
from Fig. 4 or by substitution of y into Eq. (18). Maximum on-axis gain for the 
fax field is found by adding the contributions from Fig. 11 and 12, from which we 
obtain 95.5dB and -1.5dB respectively. Referring back to Fig. 4 we can see how 
the maximum axial gain of 94 dB decreases with deviations in a from 1.07. Using 
Fig. 13 we can anticipate a 3 dB loss in gain for a .55X/D radians pointing error. 
At the 10.6 micrometer C0 2 wavelength this corresponds to a pointing error of 
29 microradians. The full cone angle (2 5) of the beam in the far field is esti- 
mated to be 92 microradians by approximating the y = 0.2 curve in Fig. 14. 


If we were to defocus the antenna or move our observation point into the 
near field the value of /5 will no longer be zero and can be calculated using 
Eq. (10). Assuming /?' = /3 /2v = 1.46, corresponding to one of the maxima in 
Fig. 5 for y = 0.2, we obtain a transmitting efficiency factor of -13.9 dB. Adding 
the contribution from Fig. 11 yields a total on axis gain of 81.6 dB. The obser- 
vation angle corresponding to the -3 dB point relative to the peak of the near field 
gain distribution is found from Fig. 8b to be ka Sin 6 - 2.5 yielding 0 - 42 
microradians. 


4.0 RECEIVER GAIN 

4.1 Loss Due to the Central Obscuration 

We now consider the gain of the previously described telescopes when used 
as receivers. We assume that the source of the radiation is sufficiently far 
away so that a uniform plane wave impinges on the receiving aperture. The 
gain of the receiver is defined as 
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where A =77 a 2 is the area of the primary aperture as defined in Eq, (11). The 
quantity (1 - 7 2 ) can be considered a receiver efficiency factor gp ( 7 ) which takes 
into account the fraction of the impinging light scattered by the central obscuration. 
The receiver efficiency factor in dB is plotted in Fig. 18. 


CM 



Figure 18. Receiver loss due to the centra! obscuration in dB. 


4.2 Loss at the Detector 

Although we have used the word ’telescope" throughout the present article, 
the primary and secondary mirrors are actually equivalent to a large aperture 
lens of focal length f. The detector is placed in the back focal plane of this 
"lens" and the aperture plane containing the central obscuration lies a distance 
d (corresponding to the separation between the primary and secondary mirrors) 
in front of the "lens". Goodman 11 has shown that the electric field in the back 
jcal plane is given by 
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where x Q and y Q and x f and y f are the Cartesian coordinates in the object plane 
and back focal plane respectively, E q is the amplitude of the uniform plane wave 
impinging on the aperture (object plane), and t () (x 0 , y (l ) is unity over unobscured 
portions of the aperture and zero over obscured areas. Converting to cylindrical 
coordinates and assuming centrally obscured apertures as before, we obtain 
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Since the '‘object 11 is circularly symmetric, the image in the focal plane is inde- 
pendent of v/ t and we can use arguments identical to those used in going from 
Eq. (3) to (4) to obtain 


eo,) - i 7' <-»i> [i j; (i -fj f] f -„'i „ 

Jl- 


k, n' < 


v r. u ' d\ n r / k ^- f 

,- f - ° XP L J 2 f V 1 * fj *! Mi- 


ka.. \ 


}\ 


i (dot 


/ * \ f / j 


where •• - b/a as before. 
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Equation (25) gives the distribution of electric field across the detector. 
We note from the exponential in r-'q, (25) that if the entrance aperture of the 
telescope does not lie in the front focal plane of our equivalent ens" the:; the 
phase front across the detector is not planar. Comparing Eq. (25; with the 
expressions in Appendix A indicates that the radius of curvature of the phase 
front is given by 
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4.2.1 Direct Detection 

In a direct detection system, the overall power into the detector is the 
important quantity. The fraction detected is simply given by 
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where R D is the radius of the detector (assumed to be circular). By defining 
the variable u - ka.^ ( /i, the above equation simplifies to 



where the fraction detected is a function of kaR and . . The fractional power 
detected is plotted in dB in Fig. 19 against the parameter kaR^/f for several 
values of . . 
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Figure 19. Froctional power detected as a function of *ne antenna obscuration ratio •yn>>d the 
quantity ka F-p/f where k = 2tt/\, a is the aperture radius of the antenna, Rp is the detector 
radius and f is the focal length of the antenna. (Direct detection or “matched" heterodyne de- 


tection described in Section 4.2.2). 


4.2.2 Heterodyne Detection 

In a heterodyne detection system, the electric field is mixed with the local 
oscillator field which, for the moment, is assumed to be a uniform plane wave 
across the small aperture of the detector. The instantaneous power into the 
detector is given, by 
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where is the real amplitude of the local oscillator field, &j t and w LQ are 
the transmitter and local oscillator frequencies, and 6 is an arbitrary lag in 
the longitudinal phase of the two beams. Substituting Eq. (25) into (29) and 
ignoring components associated with frequencies of order 2c ^ (which do not lie 
within the frequency response of the mixer), we obtain 


P(t) = 27T 


*X0 *4x> 


J (Ji(>j) -7jj(yu)] du sin [(w T - ^ L0 ) t + <*>| ( 30 ) 


where we nave also ignored the relatively minor contribution of the transmitter 
power to the D.C. component and the curvature of the transmitter wavefront at 
the detector. The curvature can be ignored if 
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which is often the case. Examination of the integral in Eq. (30) reveals that 

the detector radius R opt which maximizes the signal at the intermediate frequency 

(IF) is determined by the first root of the equation 


-yjf r 


which, as can be seen from Eq. (25), is also the determining equation for the 
nulls in the electric field distribution in the back focal plane of our equivalent 
lens. In other words, for optimum performance, the detector perimeter should 
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coincide with the first dark ring of the signal diffraction pattern. The use of a 
larger detector results in partial cancellation of the detected signal due to the 
cha T '' p e in sign of the electric field in spite of the fact that the detector captures 
mere of the transmitter energy. The general solution to Eq. (32) was obtained 
numerically and is plotted as a function of y in Fig. 20. 



Figure 20. The introduction of o central obscuration into the aperture modifies the Airy pattern 
in the plane of the detector. This curve is used to match the antenna F-number (f/2a) to the de- 
radius (R opf ) so that the detector captures only the central disk of the signal diffraction 
pattern- 


If we define P R = Eg 7ra 2 (l - y 2 )/2 as the time averaged power collected 
by the aperture and P LQ = v Rg pt /2 as dc local oscillator power absorbed 
by the detector (uniform illumination only), then the mean square value of 
power in the IF component can be written 
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( 33 ) 


The signal to noise ratio at the detector is defined as 





(34) 


where i 2 is_the mean square current generated in the detector due to the IF 
signal and^ i 2 is the current generated by a particular noise process. The 
quantity i 2 is in turn given by 



(35) 


where e is the electronic charge, t is the detector quantum efficiency, G is 
the detector gain (if any), h is Planck's constant, and v is the photon frequency. 
In most references dealing with heterodyne detection, the mean square signal 
current is generally written in the form 
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Comparing Eqs. (33), (35), and (36), we infer 



»VN VH' fW* nf 13* •?»'"v?^r7r' ffMW** > ‘i* ,,l i'* *,/»*•«.«*■,♦ v»* - 





Thus, the time average power P R actually received by the aperture must be 
multiplied by the factor ^ (y) to obtain the so-called "signal power" P s which 
appears in the standard photodetector equations. The dB value of \{y) for 
uniform illumination is plotted versus y in Fig. 21(b). The above 
calculations were performed under the assumption that the local oscillator beam 
illuminates the detector uniformly. The efficiency of the receiver can be 
improved, however, by focusing the local oscillator energy onto the detector. 

The most general system one can imagine is one in which the local oscillator 
beam is focused by a separate set of optics into the detector. If we allow the 
effective F numbers of the signal and local oscillator optics to be given by F s 
and F.^ and the respective central obscurations to be given by y s and y LQ , we 
obtain the general result 
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The latter equation is obtained by substituting two electric field distributions of 
the form given by Eq. (25) into Eq. (29) and proceeding with the calculations as 
in the uniformly illuminated case. 


A system of practical interest is one in which the detector captures only 
the central lobe of the signal Airy disk but the local oscillator beam is partially 
focused onto the detector by an unobscured lens ( 7 ^ = 0). Eq. (39) then reduces 
to 
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where kaR opt /f is defined as before and has been plotted against 7 in Fig. 20. 
Eq. (40) is plotted versus the parameter F S /F L0 in Fig. 21(a) for several values 
of 7. The limiting value of tj h as F s /F lo - 0 in the latter figure corresponds 
to unfocused uniform illumination by the LO. The comparison of these limiting 
values with the uniform illumination curve in Fig. 21(b) calculated using Eq. (38) 
at the points 7 = 0.0, 0.1, 0.2, 0.3 and 0.4 verifies this relationship. As the LO 
beam is focused, the receiver efficiency increases until the first dark rings of 
the local oscillator and signal diffraction patterns coincide. Focusing the LO 
beam still further causes the first outer intensity ring of the diffraction pattern 
to fall on the detector resulting in a partial cancellation of the signal at the 
intermediate frequency and reducing the overall efficiency of the receiver. 


The ratio of the F numbers which results in maximum efficiency for a 
heterodyne receiver is given by 



(41) 


where kaR opt /f is a function of the obscuration ratio 7 and is plotted in Fig. 20. 
The overlap integral of the two electric field distributions is maximized when the 
above condition holds and thus the signal at the intermediate frequency is strongest. 
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F$/Flo 


Figure 21g. The heterodyne detection loss in dB os a function of Fj/F^q where F $ is the 
F-number of the optical antenna and F lO is the F-number of the local oscillator focusing optics. 
The limiting values os F s /F l q Ogive the detection efficiency for uniform illumination by the 
local oscillator. The efficiency improves os the local oscillator is focused more effectively on 




the detector until the second bright ring of the LO diffraction pattern begins to fall on the de- 
tector. This results in partial cancellation of tho I F signal and reduces the detection efficiency. 
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OBSCURATION RATIO, 7 


Figure 21b, The heterodyne detection loss in dB as a function of obscuration ratio y. The lower 
boundary of the shaded region is the detection efficiency for uniform illumination by the LO. The 
upper boundary is the optimum detection efficiency obtained when the first dark rings of the sig- 
nal and LO diffraction patterns coincide with the perimeter of the circular detector. This upper 
boundary corresponds to the peaks in Fig, 20(a). The shaded region designates the intermediate 
range of efficiencies obtainable with partially focused local oscillator beams. 
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Of course, in deriving Eq. (38), we have assumed that the total DC local oscillator 
power into the detector (and hence the detector shot noise) is held constant. The 
increase in receiver efficiency brought about by focusing the LO occurs because 
we have concentrated a greater fraction of the total LO power near the center 
of the detector where the signal intensity is strongest thereby maximizing the 
overlap integral. Upon substitution of Eq. (41) into (40) we obtain 



Again, since kaR /f is uniquely defined for each value of y> the latter expres- 
sion can bo plotted versus the single parameter y to yield the upper curve in 
Fig. 21(b) labeled "Optimally Focused Local Oscillator". The latter curve 
represents a lower limit on the detector loss under the assumption that the 
detector captures onlv the central lobe of the signal distribution. The losses for 
partially focused local oscillator beams (the full central disk of the LO pattern 
does not fall on the detector) fall in the shaded region between the two curves. 

We now consider one final very special case where the F number of the local 
oscillator beam matches the F number of the signal beam (F s /F L0 = 1) an a 
mask is constructed over the focusing lens for the local oscillator such that 
y = y c = y. We obtain from Eq. (39) 

7 LO 7 S ' 
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du 


(43) 


Under the latter conditions, the electric field distributions of the two beams 
are identical in the plane of the detector except for a constant amplitude 
factor. The resulting equation for the efficiency is identical to Eq. (28) ob- 
tained for direct detection which was plotted previously in Fig. 19. Unlike the 
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situations disc issed previously, the detector dimensions can be enlarged beyond 
the central lobe in the signal diffraction pattern without encountering field 
cancellation effects . 


4.3 Receiver Summary 

Figure 11 gives the dB gain of an unobscured antenna of radius a for a 
wavelength \ . Figure 18 gives the dB loss in received power due to blockage 
of the incoming light by the central obscuration. Figures 19 and 21 give the 
dB loss at the detector for direct detection and heterodyne detection respectively. 

For a direct detection system, the total receiver gain can be obtained by 
summing the separate contributions from Figs. 11, 18 and 19. Of course, 
the effect of lossy optical elements or non-unity detector gains and quantum 
efficiencies must be separately accounted for. 

In heterodyne (or homodyne) detection systems, the effective receiver gain 
depends on the electric field distributions of the local oscillator and signal beams 
across the detector. Good performance is obtained when the first dark rings 
surrounding the central Airy disks of the focused signal and LO diffraction pat- 
terns coincide with the detector perimeter. The system experiences significant 
degradation w.V-n tho LO uniformly illuminates the detector or is only partially 
focused on the detector. Large obscuration ratios degrade the receiver gain for 
two reasons: (1) greater blockage of the incoming radiation results; and (2) a 
smaller fraction of the total received energy is contained within the central lobe. 
The total receiver gain is obtained by summing the separate contributions from 
Figs. 11, 18 and 21. The optimum ratio of detector size to antenna F number can 
be inferred from Fig. 20. The efficiency of a heterodyne detection system can be 
further improved by constructing an appropriate mask over the lens which focuses 
the local oscillator beam onto the detector. If the obscuration ratio of the mask 
matches the value of the obscuration ratio for the antenna and if the F numbers 
of the signal and LO optics are the same, then the field distributions of the two 
beams across the face of the detector are identical except for a constant ampli- 
tude factor. As a result, the detector size can be enlarged to contain the outer 
rings of the diffraction distributions without encountering field cancellation effects. 
In this very special case, the losses are identical to those suffered in a direct 
detection scheme and the total receiver gain is calculated using Figs. 11, 18 and 
19. Systems in which the local oscillator is focused onto the detector require 
more stringent alignment tolerances 12 , however, and care must be taken to insure 
that the phase fronts of the two beams do not vary widely over the detector di- 
mensions. Calculations of the power distributions in the near field and in the 
vicinity of the focus for Gaussian beams focused through annular apertures have 
been carried out by Holmes et al. 13 The present work assumes that the LO spot 
size is large compared to the effective aperture of the focusing lens. 
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4.4 An Illustrative Example 

We consider a Cassegrainian telescope with a 20-cm diameter primary and 
a 4-cm secondary mirror. Thus, we have a = 10 cm and 7 = 0.2. From Fig. 11, 
we obtain a value of 95.5 dB at a wavelength of 10.6 microns. This represents 
the theoretical upper limit on the gain available from an unobscured antenna 
having a diameter of 20 cm. Fig. 18 indicates that 0.2 dB is lost due to blockage 
of the incoming radiation by the secondary mirror. If the local oscillator 
illuminates the detector uniformly, an additional loss of -3 dB is evident from 
Fig. 21(b) yielding a total receiver gain of 92.3 dB. Focusing the LO so that 
the two central disks are barely contained within the detector raises the receiver 
gain to 94.1 dB. The appropriate ratio of detector radius to antenna focal length 
can be obtained using the relation ka R opt /f = 3.67 from Fig. 20. 
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APPENDIX A 


MATHEMATICAL REPRESENTATION OF SPHERICAL 
WAVEFRONTS AT AN APERTURE 


Consider a spherical wave emanating from a point source a distance R from 
c. diffracting aperture as in Fig. A. Consider also a Gaussian amplitude function 
superimposed on the wavefront. The electric field al the aperture can then be 
expressed as 


E 0 (r 0 ) = Ee' kf e' r o A ' 2 

Assuming the curvature of the phase front is much larger than the aperture 
allows us to write the approximation 


r 


(R 2 + r 2 ) 1 2 = R 1 + 


/ f 2 \* 2 

r r 2 1 

( 1 + — ) e- R 

1 1 

V R 2 / 

- 2R 2 _ 


The field can now be approximated by 


E 0 (r 0 ) = E o jkR e' k '* /2 * t 


-W 


where e )kK is a constant phase factor that is lost in a calculation of the intensity 
at a point beyond the aperture. 



Figure Al. Mathematical representotion of curved wavefronts at an aperture. 
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APPENDIX B 


SECOND ORDER PERTURBATION ANALYSIS OF 
OBSCURED APERTURES 


We are interested in obtaining an approximate analytical solution to the 
transcendental equation 


2q2 Z. 1 e " a2 ( 1_y2 > - i.. 

2a 2 y 2 + 1 


In Section 4, we obtained the solution a Q = 1.12 for the special case y = 0. 
We will assume that an approximate solution of the form 


a ~ a. + Sa, + 8 2 a. 


can be obtained where a. and a are the first and second order corrections to 
the zeroth order solution a Q for small y , and 8 is a dimensionless parameter 
whicj simply monitors the relative magnitudes of the individual terms. Since y 2 is 
the perturbing quantity in equation (17), we associate a factor of S with each 
factor of y 2 . Ignoring terms of order S 3 or higher allows us to rewrite (17) 

(after much algebra) as 


(2a 2 + 1) + S(4a fl a 1 ) + S 2 (4a 0 a 2 + 2a 2 ) 

= e a o jl + S(a 2 ')' 2 + 2a Q a 1 ) 

(2a a - a*y 2 ) 2 

+ £2 2a 0 a j y 2 + 2a 0 a 2 + °f + + 2 a 2 7 2 ( 2 a 0 a j _ a 2 y 2 ) 


Collecting terms of zero order in 8 gives 


(2a 2 + 1) e 0 = 1 


which is the defining equation for a 0 with the solution a Q =1.12. 
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Collecting terms of first order in 5 and solving for gives 


/ = -1.30 /. 


2(2 - e a o ) 


The terms of second order in S result in an equation for a , that is 


.2 1 2 

2a 0 (2 - e#)l 


foa^y 2 + 2 a 2 + 4a 2 a 2 + 4 / - 3a 4 /] - 


= 2.12 /. 

Thus, to second order the optimum aperture to beamwidth ratio for maximum 
axial gain for an arbitrary obscuration ratio is given by 

a-a 0 +O.J +a 2 = 1.12 - 1.30/ + 2.12/. (18) 

Comparisons of the above approximation to exact solutions of equation (13) 
obtained on a computer for several test values of y indicate an accuracy to 
within ±1% for y < 0.4. 


APPENDIX C 
PROGRAM P - GAIN 


The program, which has been named P-Gain, for power gain, is written in 
Fortran IV in the complex double-precision mode and computes near and far 
field antenna gain patterns. Eq. (12), 


&!•(&. £. y, x) = 2a 2 



-a 2 u 


J 0 (X y u) duj 


( 12 ) 


reproduced here from Section 3.1, is the basis for the calculation. The program 
owes its bulk, somewhat, to the injection of a large number of auxiliary comments 
and to the rather sophisticated plotting technique (Wolfplot) empJ •'yed. 

Care is taken to insure accuracy while attempting to incorporate speed into 
the computation. The integration is performed by a slightly modified version of 
a double precision 12-point Gaussian quadrature subroutine. The routine was 
revised to allow complex integration thereby reducing the computation time of 
the program by approximately one-half. 

The separation ot the annular region of the antenna aperture into integrable 
divisions to obtain sufficient accuracy was determined by analyzing the integrand. 
Location of the observation point, i.e. its perpendicular distance to the optical 
axis in the near or far field, dictates whether the Bessel function or the complex 
exponential function has the higher frequency of oscillation. The number of 
half cycles of the faster oscillating function prescribes the number of regions 
of integration in the program. The integration limits for each quadrature 
integration are then taken over each half cycle of that function. 
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//E 7B JKPGl JOB ( £ 7002 1 1 BHD ? T , EOOU53 »003002 )t!29ftoSGLEVEL“l P-&AIN 
it bXEC TVwTKANH 
//SOURCE. SYS1N UD * 

C THIS kROGKAM HAS the capability to calculate aNC plot POWER GAI.N CURVES 
C EUR ANALYTICALLY EXPRESSABLE CIRCULAKLY SYMMETRIC SUUKCt FUNCTION INCIDENT 
C UN ObSTRUCTEU AND UNUbSTKUCTED CIRCULAR APERTURES. THE CURVES ARE kELATlVE 
C TU A GAIN FACTOR ( **P I*A/WAVELENGTH**2 ) , WHERE A IS THE APEkToRE ARfcA. 

C THE OBSTRUCTIONS ARE ASSUME CIRCULAR AND CENTRALLY LOCATED WITHIN THE 
C APERTURE. THE INCIDENT FIELD IS ASSUMED NURMAL TO ThE APEkIURE WITH THE 

c capability of being separated into amplitude and phase components. 

C P-GAIN WAS WRITTEN SPECIFICALLY *UK LASER COmHUNC I AT I u.\S FUR THE ANALYSIS 
C OF CASSEGRA INI AN ANTENNAS. TO OBTAIN THE KtUUlKE D ACCURACY IT WAS WRITTEN 
C IN DOUBLE PRECISION. 

C all CALCULATIONS hEkE 1 -iAUE In SPHERICAL COORDINATES AS DISPLAYED IN 
C X-UUCUMENT X-52*-73~ln5. 

C 

C PROGRAM P-GAIN - BEGIN MAIN 
0 

IMPLICIT REAL*rt (B-H,0-Z) 

C THE LETTER A WILL NUT BE DEFINEO TO BE DO JHLc PRECISION hUT USED FOR SINGLE 
C PRECISION ARRAY NAMES BECAUSE WOLFPLOT PLOTTING ROUTINES WILL .,L«T ACCEPT 
C DOUBLE PRECISION ARRAYS. 

C DEFINE AS COMPLEX THE FOLLOWING - 

COMPLEX- 16 SUhEXP , FE XP , SEXP ,DCON JG 

C FEXP - A FUNCTION EXTERNAL TO THE MAIN PROGRAM AND INTcGkA»\D uF THE COMPLEX 
C INTEGRATION SUBROUTINE CD0G1Z 

C SEXP - THE RESULTANT OF THE COMPLEX INTEGRATION SUBROUTINE. 

C SUMEXP - THE SUM OF THE RESULTANTS FKUK THE CUMPLnX INTEGRATION SUbRUUT I . 
C DCOM JG - A DOUBLE PRECISION EXTERNAL FUNCTIUN CO^ JUG AT I V IG SUMEXP • 

C DIMENSION AND DEFINE AS SINGLE PRECIS1UN THE AKkaY NAMES; 

REAL** AGAIN! 20 I ) , AGA I ND ( 20 1 ) * AX I (201) 

C AX 1 ~ A PLOTTING PARAMETER FUR OFF AXIS DISTANCE WHICH ACCOUNTS FOR THE 
C beam spread. 

C AGAIN = GAIN 
C AGAIND = GAIN IN DECIBELS 

C DEFINE AS SINGLE PRECISION THE THREE VARIABLES YMlN, YSAVE, AinO XsAVE WHICK 
C SERVE AS VALUE HULDEKS DURING THE PLOTTING ROUTINE. 

REAL** YMIN,YSAVE,XSAVE 
DIMENSION NOCHAR ( 2 ) 

LOG I CAL*1 BCHAK( 8) 

DATA NOCHAR/+4,0/ 

EXTERNAL FEXP,DCONJG 
COMMON DELTA, BETA, XI, GAMMA 

c initialize the wolfplut plotting routine 

CALL PLOT ST (2001,*) 


if 


I 



3 s " 





t 


55 


"X 

s? 

A 



$ 


2f 


* ■*,! 


10 FORMAT (3015.5) 

20 FORMA f ( •O 1 »6X, ‘DELTA* f 11X, •B6TA2 1 1 10X, »X1M* X' , 10X, 1 GAMMA 1 ) 

30 FORfAT (' THE PRUGRAK HAS BEEN COMPLETED') 

40 FORMAT ( 1 0 1 * 1 X , 1 OFF AXIS DISTANCE 1 , 8X , 'GAIN* , 14X, ‘GAIN IN DB ' ) 

50 FORMAT (F13»4,6X,F13.4,9X,FI3.4 ) 

60 FORMAT (2 12,016.5) 

70 FORMAT (4015.5) 
rtT EST = 0 

C MT £ ST IS USED AS A TEST FOH THE NUMBER OF DATA SETS (MCARDS) READ 

HO HEAD ( 5*60) NCARDS, HOARDS, X1MAX 

C HOARDS = T HE NUMBER OF SETS OF DATA CARDS DICTATING THE NUMBER OF GRAPHS. 

C NCAKDS = THE NUMBER OF CARDS IN A SET OF DATA CARDS DICTATING THE NUMBER OF 

C curves per Graph. 

C X 1 MAX = K*a*SINJ THETAMAX ) WHERE A IS NOW ThE APERTURE RADIUS, K IS THE 
C PROPAGATION CONSTANT, ThETAI-.AX IS A MAXIMUM ANGLE WHICH ACCOUNTS FOR THE 
C PATTERN SPREAD, AND XlMAX IS CHOSEN AS THE ?;AX1MUM DISTANCE FROM THE AXIS 
C TH<uT ONE WISHES TO OBSERVE Tht PATTERN AND IS THE MAGNITUDE OF THE MAXIMUM 
C VALUE OF THE ABSCISSA. 

A X If* AX = X 1MAX 
AX1MIN=-X1MAX 

C AX 1 MA X AND AXlMIN ARE MAXIMUM AND MINIMUM X-AXIS PLOTTING PARAMETERS. 
OT£ST=0 

C •'“TEST IS USED AS A TEST FUK A SET OF DiTA CARDS (NCAKDS) TO BE READ IN 
C IhUI Vli^ALL*- EACH DATA CAkD CONTAINS VALUES FOR DELTA, BETA, AND GAMM*. 
hTEST =MTEST + 1 
CALL CuNDNS 

C SET UP SUBJECT SPACE FUk PLOTTING 

CALL OGRID ( AX 1M I N , AX 1 MAX , 10 , • F d. 1 * , 1 . -40.0,0. 0, H, 1 F5 • 1 »,1,U) 

C COMPLETE TICK PARKS ALONG Tub TUP BORDER OF THt GrvAPH 

CALL OGRID ( A X 1 M I N , AX 1 MAX , 10 , • F 5. 1 1 ,0, 0.0, -40. 0,0, 1 F5. 1 *«0,l0) 

C COMPLETE TICK MARKS ALONG T H t RIGHT BORDER UF THE GRAPH 

CALL OGRID ( AX1MAX, AXIMN,0, *F5.1 ' , 0 , -40 . 0 , 0. 0, 8 , 1 F 5 . 1 *,0,10) 

90 READ (5,10) DtLTA,rt£TA2 * GAMMA 

C DELTA = A/NO, WHERE A = APtKTURE RADIUS AND WO = THE GAUSSIAN BEAM WIDTH. 

C bFTA2 = BETA/2PI 

C BETA = (K*A**2/2)*( 1/Rl+l/K), WHERE K = 2*P I /WAVtLEiNGTH , R=THE DISTANCE 
C HETWEEN A POINT IN THE APERTURE AND THE OBSERVATION POINT, AND K1 IS THE 
C DISTANCE FROM THE SOURCE TO THE POINT IN THE APERTURE. 

C GAMMA = 8/A - THE FRACTION OF the APERTURE TO BE OBSCURED, WHERE b = THE 
C RADIUS OF i HE CENTRAL OBSCL-' A T 1 f N. 

NT E SI GTFST + 1 
WRITE (6,20) 

WRITE (6,70) DELI A , MET A2 fXi,.AX* GAMMA 
WRITE (6,40) 

DO 1 7( 1=1,101 





FSUB=I— 1 

C DEFINE THE X-AXIS ARRAY AX1 

XI=X1MAX-FSUB*X1MAX/100.D0 

AXU I )=X1 

K=202- I 

AX1 ( K )=— XI 

PI=3. 141592600 

BETA=BETA2#2*p I 

C CALCULATE The INTEGRAL NUMBER OF HALF CYCLES IN THE INTERVAL (0,1) BY (1) 

C. DIVIDING THE ARGUMENT OF THE COMPLEX EXPONENTIAL BY PI, WHICH IS Nl, OR (2) 

C DIVIDING THE ARGUMENT OF THE BESSEL FUNCTION BY PI, WHICn IS N2. 
N1=BETA/PI+1.D0 
N2=X 1/PI+ 1 , DO 

C WHICH OSCILLATES FASTER I THE EXPONENTIAL OR THE BESSEL FUNCTION) WILL 
C DETERMINE THE NUMBER OF CYCLES OVER WHICH TO CARRY UUT THE INTEGRATION TO 
C ASSURE ACCURACY OF THE SECONDARY PATTERN. FARTHER FROM THE AXIS THE BESSEL 
C FUNCTION IS ANTICIPATED TO OSCILLATE FASTER AND WILL THEREFORE RESULT IN A 
C LARGER NUMBER OF HALF CYCLES TAKEN WITHIN THE ANNULUS OF THE APERTURE. 
NHCYC=KAX0(Nl,il2) 

C INITIALIZE THE LUWfcR AND UPPER INTEGRATION LIMITS RL AND RH, RESPECTIVELY. 
RL=GAHMA<-*2 

IF (NHCYC.E0.N1) GO TU 100 

C IF XI EQUALS ZERO (THOUGH NOT LIKELY) AN ERROR WILL BE DETECTED IN COMPUTING 
C RH IN THE NEXT STEP, THEREFORE TU AVOID THIS A TEST IS MADE FIRST. 

IF (Xl.EO.O.DO) GO TO 110 

RH=(PI/X1)**2 

GO TO 120 

C IF BETA EQUALS ZERO AN ERROR WILL BE DETECTED IN COMPUTING KH IN THE NEXT 
C STEP, .HEREFURE TO AVOID THIS A TEST IS MADE FIRST. 

100 IF ( BETA. EU. 0,00)' GO TO 110 
RH=PI/BETA 
GO TO 120 
110 RH=1 , DO 
120 CONTINUE 

C INITIALIZE INTEGRAL VALUES TO ZERO 
SUMEXP= (0.00,0.00 
C COMPUTE THE INTEGRALS 

C THE INTEGRATION IS CARRIED OUT OVER EACH HALF CYCLE ANU SUMMED TOGETHER 
C IN SUMEXP IN THE FULLOW IMG DU LOOP. 

DO 150 J=1,N8CYC 

C TEST TO SEE IF LOWER LIMIT IS GREATER THAN UPPER LIMIT. IF IT IS, THE 

c program jumps over a number of steps to readjust rh. 

IF (RH.LE.RL) GO TO 130 
IF (RH.GT.l.DO) RH* 1 . DO 
CALL CDQG12 ( RL , RH, FEXP , SEXP ) 





! 

1 

) 



SUMEXP=$UMEXP+SEXP 

C SET LOWER LIMIT EQUAL TO THE UPPER LIMIT AND ADD ONE MORE HALF CYCLE TO RH. 
RL=RH 

130 IF (NHCYC.E0.N1) CO TO 140 

C CABLED LINE 23 TESTS FOR THE METHOO by WHICH NHCYC WAS CHOSEN (N1»N2) so 
C AS TO OIRECT THE PROGRAM TO INCREMENT BY ONE HALF CYCLE THE CORRECT UPPER 
C LIMIT RH. 

IF (Xl.EO.O.UO) GO TO 160 
RH*!DSORT( RHI+PI/Xl )**2 
GO TO 150 

140 IF ( BETA .EO.O .DO ) GO TO 160 
RHsRH+PI /BETA 
150 CONTINUE 

C LABLED LINE 21 SEP VES AS A JUMPING POINT OUTSIDE THE ABOVE LOOP TO AVOID 
C REDUNDANT CALCULATIONS WHEN Rl.=RH» 

160 CONTINUE 

C COMPUTE THE VALUE OF THE GAIN 

8GAlN=2.D0*(DELTA**2)*SUMfcXP*0C0NJG(SUMEXP) 

AGAIN! I ) =BGA IN 
AGAIN! K)=BGAIN 
170 CONTINUE 

C CONVERT THE GAIN DISTRIBUTION TO 0EC1BELS AND DEFINE THE Y-AXIS ARRAY AGAIND 

00 1B0 1=1.201 

AGAIND! I )=10.*AL0G10( AGAIN! I ) ) 

180 WRITE (6,50) AX1 ( I ) .AGAIN! 1 ) .AGAIND(I) 

C The PROGRAM FROM THIS POINT TO LABEL 17 PLOTS THE GAIN DISTRIBUTION IN DB 
C AND EXTRAPOLATES THE CURVE TO YM1N WHEN AGAIND FALLS BELOW YMIN. 

YMIN»-40.0 

1 = 1 

IF ( AGAIND! 1) .LT.YHIN) GO TO 210 
C PLOT GAIN DISTRIBUTION ABOVE YMIN 
190 10=1 

CALL INCORE ! BET A2 .BCHAR , 15, 1 ,5, 2,0 ) 

00 200 1=1,200 

IF (AGAIND! 1+1 ).Gfc. YMIN) GO TO 200 
XSAVE = AX1 ( 1 + 1 ) 

YSAVE = AGAINl>! 1 + 1) 

AXl! 1 + 1 ' = AX1( I ) + ( AXl! I + D-AX1! ! > )*< YMIN-AGAi;.U( 1 ) )/ l AGAIND! 1 + 1 l-AG 
1AIND! I ) ) 

AGAIND! 1+1 )=YMIN 
AHGLD=75.-10 

CALL PLOTS (AXl! 10) .AGAIND! 10 ), 1+2-10,0.0,0, 0.0,0, 0.0,0, l ,0.0,0. 0, 

10.0,0 .BCHAR! 2) , NOCHAR, AHOLD, 0,0) 

1=1 + 1 

AXl! I )=XSAVE 








AGA 1 NO ( I)=YSAVE 
GO TO 210 
200 CONTINUE 

C THIS PLOTS THE' TAIL OF THE GAIN DISTRIBUTION ABOVE YHIN OK THE WHOLE 
C DISTRIBUTION IF AGAINO NEVER DROPS BELOW YM IN* 

AHOLD=70.-I0 

CALL PLOTS (AX1( lOI.AGAINDI 10 ) , 202- 10,0.0,0, 0.0, 0* 0.0, 0, 1 ,0.0, 0.0, 
10.0,0,l«CHAK( 2) , NOCHAR, AHULI', 0,0) 

GO TO 230 

210 DO 220 1*1,200 

IF (AGAlNOI 1*1) .LT. YHIN) GO TO 220 

AX1 1 I )=AX1! 1 >♦( AX1< I + U-AXH I » )*< YMIN-AGAINOI 1 )) / ( AGA INO( 1 + 1 I-AGAI 
1ND( 1 ) t 

AGAINOI 1 >=YHIN 
Gil TO 190' 

220 CUNTINUE 

C RETURN FOR NEW DATA CARD 

230 IF (NTuST.LT.NCAKDS) GO TO 90 

C THE NEXT FOUR CALL STATEMENTS ARE TO LABEL ThE AXES AND ADVANCE THE FRAME 
C TO PLOT A NEW SET OF CURVES. 

CALL ROTLET 11.57081 

CALL VERLIN l')B0( RDTCAF N 1 AG REWOP » , 22 , . 25 * 5. 5,0, 0) 

CALL ROTLET (0.0) 

CALL H0RL1N PX - K*A*S IN ( ThET A ) 1 ,19, 5. 25,. 60, 0,0) 

CALL FKMADV 

C RETURN FOR NEW DATA SET 

IF ( MTEST .LT.MC ARDS ) GO TO 80 
CALL SETHLK (999) 

CALL ENDPLT 
WRITE (6,30) 

RETURN 

FNO 

C COMPLEX DOUBLE PRECISION GAUSSIAN INTEGRATION SUBROUTINE. 

SUBROUTINE CDC1G12 ( XL »XU,FCT ,Y ) 

IMPLICIT REAL*B (A-H,0-Z) 

C0MPLEX*16 FCT , Y 
A=.500*( XU+XL ) 

H=XU-XL 

C=.4907803171233596?D0*B 

Y=. 2.3587668 1932 559 140-1*1 FCT ( A+C >+FCT( A-C) ) 

C=. A 52058628 1852374300*B 

Y =Y+. 5346966299765921 5D-1*< FCT ( A+C ) +FCT ( A-C ) ) 

C = . 3 8495 133709 71 523400*8 

Y=Y+. 80039 16427 16731 10-1* ( FCT ( A+C )+FCT< A-C )) 

C=. 29365 89771433087200*8 


S 

i\ ' 

I 


*5 
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Y=Y+ .101 583713361 5329600* (FCT( A+C>+FCT( A-C ) ) 

C=. 1839 15 7494990901 000* B 

Y=Y+. 11674626826917740D0*( FCT( A+C )+FCT ( A-C ) ) 

C=. 6261 67042557344580-1*8 

Y = B* ( Y+. i. 24573 52 2 906 701 3900* t FCT ( A+C ) +FCT ( A-C ) ) > 

RETURN 

END 

C DOUBLE PRECISION BESSEL FUNCTION SUBROUTINE. 
SUBROUTINE DBESJ ( X ,N ,B J ,0 * I ER ) 

IMPLICIT REAL*8 <A-H,0-Z) 

B J=0. DO 

IF (Nj 10,20,20 
10 IER-1 
RETURN 

20 IF (X) 30 ,30,40 
30 I ER=2 
RETURN 

40 IF (X-15.D0) 50, 50*60 
50 NTE$T=20. 00+10. D0*X-X**'2./3. ) 

GO TO 70 

60 NTEST=90.D0+X/2.D0 
70 IF (N-NTEST) 90,80,80 
80 188=4 
RETURN 
90 I ER=0 
N 1 =N + 1 
BPRE V=0 .00 
C 

C COHPUTE STARTING VALUE OF M 

C 

IF (X-5.D0) 100,110,110 
100 MA=X+6.00 
GO TO 120 

110 MA=1.4D0*X+60.D0/X 
120 I X=X 

MB=N+ ! X/4+2 
MZEKO*KAXO(MA,hB) 

C 

C SCT UPPER LIMIT OF M 

C 

MRAX=NTEST 

OU 210 M»MZ6KU,MhAX,3 
C 

C SET F ( M ) , F ( M- 1 ) 

C 




h 


- 4 , 



FMl=l,00-?8 

FM=0.DO 

ALPHA=O.DO 

IF <M-IK/2>*2) 140*130 , ) AO 
130 JT=-1 

GO TO 150 
140 JT=1 
150 h2=M-2 

DO 180 K-1.M2 
MK=M-K 

BMK=2.D0*MF*FM1/X-FM 

FM*FM1 

FM1=BMK 

IF IKK-N-1) 170.160,170 
160 BJ*BMK 
170 JT*-JT 
S = 1*JT 

180 ALPHA*ALPHA+BMK*S 
BMK=2.D0*FM1/X-FM 
IF INI 200,190,200 
190 BJ=BHK 
200 ALPHA =ALPHA+BMK 
BJ*BJ/ALPHA 

IF ( G ABSC BJ-BPRE V I -DAbS ( 0*B J ) ) 220,220,210 
210 BPKEV*HJ 
IER=*3 

220 RETURN 
END 

C DOUBLE PRECISION FUNCTIUN CONJUGATING THE COMPLEX INTEGRATION RESULT ANT 
C SEXP. 

Complex function oconjg*16 iargi 

C0MPlEX*l6 ARG 
REAL*B X,Y 
X«ARG 

Y*(0«00, l.DO)*ARG 

OCONJG»X*I l. DO, 0.00I4Y*! 0.00,1.00) 

RETURN 

END 

C DEFINE THE COMPLEX FUNCTION FEXP - THE INTEGRAND UF C0UG12. 

COMPLEX FUNCTION FEXP*16 (U) 

IMPLICIT RE AL*H IA-H.0-2) 

C0MPLEX«16 ARG 
COMMON DELTA, BETA, XI, GAMMA 
10 FORMAT ( 10X, 12 ) 

X»X1*0S0RT(U) 
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EPS=1.0-05 
IF (X-EPS) 20,20.30 

20 BJ = 1 .000 
GO TO 50 

30 CALL DBESJ ( X ,0 ,B J ,EPS , I ER ) 

IF t IER.NE.O) GO TO 40 
GO TO 50 

40 WRITE (6,10) IER 

C AMP = THE AMPLITUDE OF THE INCIDENT SOURCE FUNCTION. 

50 AMP=DEXP(-(DELTA*#2)*U) 

C DEFINE THE ARGUMENT OF THE EXPOhcNTIAL AS COMPLEX. 
ARG=(0.D0,1.D0)*BETA#U 
FEXP=CDEXP(ARG)#AMP*BJ 
RETURN 
END 

/* 

// EXEC LOADER, REGION. G0=200K 

//GO.SYSLIB DO DSN=SYS2.WOLFPLOT,DISP=SHR 

//GO. PLOT APE DD DCB=( ,DEN=1 ) ,LABEL=( , BLP , , OUT ) ,UN I T= ( 7TR ACK, , DEFER ) , 

// DSN=8JKKLEIN*V0L=SER=( TO 5600 ) 


//GO. 

DAT A5 DD * 



410 

20,000 




1.121D00 

1.04000 

0.0000 


1,12 1DO0 

1.4000 

0.0000 


1.121D00 

2.02000 

0.0000 


1.121D00 

2.44000 

0.0000 

410 

20.D00 




1.121000 

0.0000 

O.ODOO 


1,12 IDOO 

0.52000 

0.0000 


1,121000 

1.22000 

O.ODOO 


1.121D00 

1.71000 

0.0000 

410 

20.000 




1.109000 

1.05000 

0.1D00 


1.109000 

1.41000 

0.1000 


1.109000 

2.04000 

0. 1000 


1.109D00 

2.47000 

0.1D00 

410 

20.000 




1.1C9D00 

O.ODCO 

0.1D00 


1.109000 

0.53000 

0. IDOO 


1.109000 

1.23000 

0.1000 


1.109000 

1.73000 

0.1000 

410 

20.000 




1.074000 

1.08000 

0.2000 


1.074000 

1.46D00 

0. 2000 


1.074D00 

2. 1000 

0.2000 



«? 





4 »» 


t . 


f . 





1*074000 

41 f > 20*000 

1*074000 
1 . 074 DO o' 

1 .074000 
1.074000 

410 20*000 

1 *028000 
1*028000 
1.028000 
1.O28DO0 

410 20.000 

1.028000 
1 . 0 2 H U 0 0 
1*028000 
1.O28D00 

410 20.000 

{). 976000 
0.976000 
0.976000 
0.976000 

410 20.000 

0.976000 
0,976000 
0 • 9 76000 
0.976000 

/# 


2.55000 

U • 2000 

0. 01)00 
0.54000 
1.27000 
1.78000 

0.21)00 

0.2U00 

0.2000 

0.2000 

1.1 3DOO 
1.55000 
2.21000 
2 . 69000 

0.3000 

0.3000 

0.3000 

0.3000 

0 .OOUO 
0.57000 
1.34000 
1 , bbDOO 

0.3000 

0.3000 

0.3000 

0.3000 

1*21000 

1.69000 

2.39000 

2.92000 

0 • 40 Ou 
0.4000 
0.4000 
0.4DOO 

0.00000 
0.6 1000 
1.45000 
2.04000 

0.4000 

0.4000 

U.4D0U 

0.4000 


S 

{. 

> 

X 







* 

* 




NASA - G 3 FC 


63 




